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iThe resistivity for an extended two channel Anderson 
model is calculated. When the temperature decreases, it 
decreases logarithmically, and has T 1 / 2 anomaly at very 
low temperatures, as seen in some dilute U alloys. The 
single-particle and the spin excitation spectra are also cal- 
culated. The low energy properties and the temperature 
dependence of the resistivity in the low temperature re- 
gion are well described by the energy scale deduced from 
the zero-energy limit of the spin excitation spectrum. 



'The two channel Kondo model (TCKM) has been stud- 
ied extensively since Nozieres and Blandin pointed out the 
possibility of the non- Fermi liquid (NFL) behaviours of its 
ground state [1]. The analytic solutions for the thermo- 
dynamic quantities have been obtained by means of the 
Bethe-Ansatz method [2]. The 7-coemcient and the sus- 
ceptibility, x, diverge as —InT at low temperatures. 

Cox has noted that a model for the electronic state of 
ions with /^configuration, such as U 4+ , can be mapped 
On the TCKM when the crystalline field ground state is 
the non-Kramers doublet [3]. Many studies have been 
reported on the applicability of the TCKM to real systems 
[4-8]. 

In experiments, a number of dilute U-ion alloys have 
been known to show the NFL behaviours. Amitsuka and 
Sakakibara have shown that 7 and x °f Thi_ x U x Ru2Si2 
have the —InT divergence in the dilute limit. At the 
same time, the resistivity decreases logarithmically with 
decreasing temperature, and then decreases in the form 
-R(O) — crT 1 / 2 with negative cr at very low tempera- 
tures [9]. 

Usually cr is expected to be positive for the TCKM. 
Though Affleck et al. have noted that cr can be negative 
in the strong exchange cases [10], recent studies based on 
the numerical renormalization group (NRG) method have 
shown that the coefficients of the logarithmic divergence 
of 7 and x are very small for such model [7,11]. Therefore 
it seems inconsistent with the experimental results. 

In a previous paper, we have presented an extended An- 
derson model which has the possibility of showing the di- 
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vergent susceptibility and decreasing resistivity with de- 
creasing temperature [12]. The single particle excitation 
spectrum pf(uS) is given as Pf(oj) = Pf(0) — Cj\io\^l 2 in 
the u> — > limit. The coefficient c/ is negative in the 
weak exchange cases and positive in the strong exchange 
cases contrasted to the TCKM. Therefore, the resistiv- 
ity is expected to decrease with decreasing temperature in 
the weak exchange cases. In this letter, we calculate the 
temperature dependence of the resistivity for the model, 
and show that it actually decreases at low temperatures. 
In addition, we show that the temperature dependence of 
the resistivity, the spin excitation spectrum, p s (uS), and 
the single particle excitation spectrum, pf{ui), are well de- 
scribed by one energy scale in the low energy region. The 
energy scale is extracted from the low energy limit of the 
spin excitation, p s (0). 

The extended Anderson model we consider is, 
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The operator f m a(ckma) is the annihilation operator of 
the localized /-electron (conduction electron with wave 
number fe). The channel is denoted by m, which takes 
and ±1, and each has spin freedom, a = ±1. The quanti- 
ties U, J and Ek are the Coulomb interaction constant, the 
exchange constant and the conduction band energy, re- 
spectively. We assume that the band has width from — D 
to D with D = 1. The density of states of the conduction 
electrons, p, is chosen as constant p = 1/2D. The energy 
level of the 0-channel, e /o is assumed to be deep enough so 
that nfo is always restricted to 1. In addition, we assume 
that the energy levels and hybridization matrices of other 



1 



two channels are independent on m, i.e. e/ TO = Sf and 
V m = V for m = ±1. If the orbits of the 0-channel and 
the ±1 channels are ascribed, respectively, to the doublet 
orbit and the quartet r 8 orbit in the cubic crystalline 
field, the present model is similar to the original Ander- 
son model of Cox, from which he derived the TCKM based 
on the Schrieffcr- Wolff transformation and restricting the 
effective manifold to /^configuration [13]. 

In this study, we calculate the electric resistivity, R(T), 
neglecting the vertex correction term. The normalized re- 
sistivity is given by 



R(T)/R(0) = 



£oi(0) 
Loi(T)' 



(7) 




where L m i is defined as 

r+oc / qj 
duo 



L m l(T) 



uj m T {uj,T) l duj. (8) 



Here, the relaxation time for the conduction electrons, 
t(oj,T), is given by using the single particle excitation 
spectrum of /-electron pf(u),T) as 



r(w,T) = 



(nV*p)p f (u;,Ty 



(9) 



The excitation spectra at finite temperatures are calcu- 
lated by using the NRG [14,15]. Parameters are chosen 
mostly to be the same as in our previous paper [12]. The 
spin excitation spectrum, p s (u>), is defined as the imag- 
inary part of the susceptibility of the operartor, Sj 
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■J m . The channel excitation spectrum, p c (u>), is 
= E ra ,iL«,±i™/™«' We define an energy 
scale, Xjf, which corresponds to the Kondo temperature 
of the fictitious two channel Anderson model without the 
0-channel spin (FTCAM)[15,16]. 

Figure 1 shows p s (v) and p c {w) for various J from the 
weak exchange cases, J/T^ < 0.35, to the strong exchange 
cases, J/TQ > 0.35. In the weak exchange cases (Figs. 1A 
and B), p s (u)) has two different structures: one at very 
low energy and another at about X^. Moreover, p s (w) has 
the a finite value in the u> — > limit, in Fig. IB. It shows 
the plateau and decreases steeply at the energy character- 
ized by 0.01/p s (0) ~ 10~ 5 as u increases. The finiteness 
of p s (0) is contrasted to the Fermi liquid behaviour go- 
ing to zero proportionally to u> in the u — > limit, and 
it indicates the — InT divergence of the spin susceptibil- 
ity [17]. As the exchange coupling, J, decreases, this NFL 
structure moves to the lower energy region and p s (0) in- 
creases (Fig. 1A). In the strong exchange cases (Fig. ID), 
p s (0) x T£ seems to be always less than 1.25 [12]. In Fig. 

IA, there is a distinct peak at uj ~ T^(= 0.19). In Fig. 

IB, the peak merges into the tail of the NFL structure, 
and looks like a shoulder. In the weak exchange cases, the 
height of this peak is much smaller than that of the NFL 
structure. When J is small enough so that the peak is 
isolated, it is quite similar to the peak in p c (oj) (Fig. 1A). 
These peaks in p s (w) and p c {oj) are formed by the Kondo 
effect of the FTCAM [16]. 



Figure 1: The spin and channel excitation at T = as 
a function of the logarithm of the energy, for various J. 
U = 1.0, s f = -1.5 and V 2 p = 0.1. The left ordinate is 
the scale for the solid line which shows p s {oj). The right 
ordinate is the scale for the dashed and dot-dashed lines 
which show p s (uj) and p c {uo) in the magnified scale. The 
Kondo temperature, T^, for the fictitious two channel An- 
derson model is estimated to be w 0.19. The sign of the 
coefficient for the M 1 / 2 term in the single particle excita- 
tion, cf, is negative for A (less than — 10 4 ) and B(— 33), 
almost zero for C (—0.3), and positive for D (0.65). 
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Figure 2: The normalized spin excitation spectra as a 
function of the normalized energy u> x p s (0), for various 
J and V. U = 1.0 and Ef = —1.5, other parameters are 
shown in the figure. For the all cases, the sign of c/ is 
negative. 



In Fig. 2, the normalized spin excitation spectra, 
Ps(u))/p s (0), are shown as a function of u>p s (0). The spec- 
tra show good overlapping in the cases of large p s (0). 
When p s (0) is large enough so that the NFL structure 
is well separated from the p c -like peak at w ~ T£, the 
NFL structure seems to be characterized by one energy 
scale deduced by p s (0). 

The single particle excitation spectra, Pf(oj), at T = 
are shown in Fig. 3. When u decreases, increases 
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Figure 3: The single particle excitation spectra as a func- 
tion of energy at T — 0, for various parameters. The main 
figure shows the spectrum in the energy range w X^- The 
inset figure shows the cusp structure at around the Fermi 
energy. The abscissa of the inset is in the square-root 
of the energy. The parameters are: solid line (U = 1.0, 
J = 0.0075, V 2 p = 0.05, T° = 0.11, p s (0) = 1.02 x 10 3 ), 
dotted line (1.0, 0.01, 0.05, 0.11, 119), dotted-dashed line 
(1.0, 0.02, 0.1, 0.19, 625), dashed line (1.0, 0.03, 0.1, 
0.19, 29.7), two-dotted-dashed line (0.1, 0.04, 0.3, 1.04, 
2.50 x 10 7 ). 
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Figure 4: The normalized single particle excitation as 
a function of the normalized energy wp s (0). The main 
figure shows p/(w), the abscissa is in the square-root of the 
energy. The straight fine line represents a fit to Pf{w) = 
Pf(0) — c/IcjI 1 / 2 . The inset figure shows same quantities 
in the logarithmic energy scale. The parameters for each 
lines are same as those in Fig. 3. 

in the energy region u> > almost analogously to the 
spectrum for the FTCAM. It shows the peaks in both 
sides of the Fermi energy, and decreases to half of the 
intensity for the FTCAM as u> goes to zero [12]. In the very 
low energy region, Pf{ui) has the cusp singularity which is 
described as Pf(uj) = p/(0) — c/|w| 1//2 . As J decreases, the 
width of the cusp structure becomes narrow. 

Figure 4 shows the normalized single particle excita- 
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Figure 5: The normalized resistivity as a function of nor- 
malized temperature Tp s (0). The parameters are: solid 
line (U = 0.1, J = 0.04, V 2 p = 0.3, p s (0) = 2.50 x 10 7 ), 
dashed line (0.5, 0.10, 0.3, 182), one-dotted-dashed line 
(0.1, 0.10, 0.3, 53.4), two-dotted-dashed line (1.0, 0.03, 
0.1, 29.7). The fine line represents a fit to R(T) = 
R(0) — crT 1 / 2 . The dotted line is the fictitious resistivity 
which is calculated without the temperature dependence 
of p f (u). The parameters are (0.1, 0.04, 0.3, 2.50 x 10 7 ), 
same as those of the solid line. 

tion, pf(oj)/pf(0) at T = 0, as a function of wp a (0). The 
low energy part seems to have a universal shape when 
we use l/p s (0) as the energy unit. In the energy region 
uj < 0.1/p s (0), Pf(u>) decreases with decreasing u almost 
linearly with |w| 1//2 . Moreover, the universality seems to 
be valid in the higher energy region where the curves de- 
viate from the simple |(jj| 1//2 law. 

The electric resistivity is shown in Fig. 5. As ex- 
pected from the numerical results of pf(u>), the resistiv- 
ity decreases with decreasing temperature. The temper- 
ature dependence of the decreasing behaviour is well de- 
scribed by using l/p s (0) as the temperature unit. How- 
ever, we note that the resistivity increases with decreas- 
ing temperature at higher temperatures T ^S> T§_, simi- 
larly to the resistivity of the usual Kondo effect. At low 
temperatures T < 0.05/p s (0), the resistivity is given as 
R(T) ^ R(0) - cbT 1 / 2 , with negative c R . Affleck et d. 
have predicted the relations between the coefficient for the 
T 1 / 2 term and that for the — InT divergent term in ther- 
modynamic quantities for the TCKM [10]. These relations 
will be checked for the present model in the future paper. 
In the middle temperature region, 0.05 < Tp s (Q) < 1, 
R(T) seems to have the logarithmic temperature depen- 
dence. 

The fictitious resistivity for which we use pf(u>,T = 0) 
instead of pf{u),T) in eq. (1) is shown in Fig. 5. Though 
it is slightly smaller than the proper resistivity, its temper- 
ature dependence is not different much from that of the 
proper resistivity [15]. 

Figure 6 compares our numerical result of the resistiv- 
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Figure 6: The comparison of the resistivity between our 
numerical result and the experimental result. The circles 
show the experimental value for Th .97Uo.o3R.u 2 Si2. The 
fitting function is R B + A(R(T) / R({))) . The parameter 
Rb is ascribed as the residual resistivity constant of the 
ThRu2Si2- The dashed line is given by fitting resistivity 
at 10K, and the dotted-dashed line is given by supposing 
2.6pf2cm of the plateau as the unitarity value plus Rb- In 
each case, l/p s (0) is determined by fitting the T 1 / 2 term 
at very low temperatures. 

ity with that of Thi_ x U x Ru 2 Si2 [9]. The calculated result 
is qualitatively similar to the experimental result. But 
the ambiguities for determination of the parameters re- 
main [18]. Further study is needed to give a quantitative 
conclusion [19]. The present model has been by no means 
directly related to the realistic model of U ion, but the cal- 
culation of the model will support the scenario to explain 
the anomalous properties based on the NFL behaviours of 
the TCKM type. 

In summary, we have calculated the dynamical excita- 
tions and the electric resistivity for an extended two chan- 
nel Anderson model. The leading temperature dependence 
of the resistivity is proportional to T 1 ' 2 at very low tem- 
peratures. Also there is the temperature region where 
the resistivity decreases logarithmically. The behaviours 
are qualitatively similar to the experimental results for 
Thi_ x U x Ru 2 Si2. We have shown that the low energy 
properties are well arranged by one energy scale deduced 
from the quantity p s (0), when the normalized exchange 
coupling is small, J/Tg < 0.35. As J/T£ increases, the 
energy region characterized by the two channel Kondo ef- 
fect spreads to the higher energy side. The sign changing 
of Cf seems to occur when the spread of the TCKM region 
becomes comparable to T^. 

In this letter, we reported only the symmetric case which 
the occupation number is n/ = 3 accounting the local 
spin. Even for the small occupation number cases, n/ < 3, 
the T dependence of the resistivity at low temperatures is 
essentially equivalent to that of the symmetric case. Detail 
will be shown in the future paper. 
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